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Abstract. We present a low energy-theory for non-linear transport in finite-size interacting single-wall 
carbon nanotubes. It is based on a microscopic model for the interacting p z electrons and successive 
bosonization. We consider weak coupling to the leads and derive equations of motion for the reduced 
density matrix. We focus on the case of large-diameter nanotubes where exchange effects can be neglected. 
In this situation the energy spectrum is highly degenerate. Due to the multiple degeneracy, diagonal as 
well as off-diagonal (coherences) elements of the density matrix contribute to the nonlinear transport. At 
low bias, a four-electron periodicity with a characteristic ratio between adjacent peaks is predicted. Our 
results are in quantitative agreement with recent experiments. 



PACS. 73.63.Fg Nanotubes 
single-electron tunnelling 



71.10.Pm Fermions in reduced dimensions - 73.23.Hk Coulomb blockade; 



1 Introduction 

Single-walled carbon nanotubes (SWNTs) are graphene 
sheets rolled up to seamless cylinders, which possess spec- 
tacular mechanical and electrical properties [1,2J. In par- 
ticular, as suggested in the seminal works |3,4,5J, due to 
the peculiar one-dimensional (ID) character of their elec- 
tronic bands, metallic SWNTs are expected to exhibit Lut- 
tinger liquid behavior at low energies, reflected in power- 
law dependence of various quantities and spin-charge sep- 
aration. Later experimental observations have provided a 
confirmation of the theory [BUT] . 

In order to study the internal electronic properties, 
including the effects of electron - electron correlations, a 
quantum dot setup, where a SWNT is coupled weakly to a 
source and drain contact as well as capacitatively to a gate 
electrode (cf. Fig. [I]), is a very well suited device. In fact, 
the current as a function of the applied bias voltage Vb 
and gate voltage V g depends on the energy spectrum but 
also on the actual form of the eigenstates of the consid- 
ered system, which for example determines whether some 
transitions are allowed or forbidden. At low enough tem- 
peratures and bias voltage, adding new particles — and 
therefore transport across the quantum dot — is hindered 
by the Coulomb repulsion [8]. Eventually the energy cost 
for placing an electron in a given energy level adds to the 
effect of Coulomb blockade. In metallic SWNTs two bands 
cross at the Fermi energy. Together with the spin degree 
of freedom this leads to the formation of electron shells, 
each accommodating up to four electrons. 

As a result, a characteristic even-odd [9j or fourfold 
[T0 l fTT | fT2] periodicity of the Coulomb diamond size as a 



function of the gate voltage is found. While the Coulomb 
blockade can be explained merely by the ground state 
properties of a SWNT, the determination of the current 
at higher bias voltages requires the inclusion of transitions 
of the system to electronic excitations. For an interacting 
ID system the occurrence of spin-charge separation is ex- 
pected, i.e. the collective excitations can be divided into 
independent modes with different spectra. In the case of 
metallic SWNTs one finds three so called neutral modes 
[4], whose energies coincide with the corresponding ones 
of the noninteracting system. Of the three modes two can 
be identified as collective spin excitations and one as col- 
lective charge excitations. Additionally, there is one mode 
of charge excitations with energies increased compared to 
the neutral modes due to the repulsive electron - electron 
interaction. In [TT|[T2] not only Coulomb diamonds but 
also low lying excitation lines have been resolved in I - Vb 
- V g measurements. Since the associated energies are not 
large enough in order to identify some of the observed 
lines with transitions to excitations of the interaction- 
dependent mode, the position of all the measured lines 
should be explainable by invoking the excitation spectrum 
of the system with neutral modes only. 

Additionally in [12] the effect of an applied magnetic 
field on the transport properties of a SWNT quantum dot 
was examined. Clear evidence for an exchange splitting of 
the ground states with two electrons in the highest occu- 
pied shell was found. Exchange effects were also inferred 
from the excitation lines of the spectrum of samples A and 
B of [11], in agreement with predictions of the mean-field 
theory presented in [13j . On the other hand, the measured 
data of sample C in [IT] do not show any sign of exchange 
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Fig. 1. SWNT quantum dot setup. A metallic SWNT is cou- 
pled via tunnelling junctions to a source and a drain contact. 
The electrochemical potential in the dot is adjusted by a gate 
voltage. 



effects, unless an unreasonably high exchange energy is as- 
sumed. The magnitude of the exchange effects depends on 
the spatial extension of the electron wave functions and 
scales like l/-/V at oms, where iV a toms is the number of car- 
bon atoms in the particular SWNT [15]. Therefore, e.g. 
the ratio between the level spacing £o of a noninteracting 
SWNT, scaling like and the exchange related ener- 
gies is roughly proportional to £>, where L and D are the 
nanotube length and diameter respectively. On that score 
it is interesting to note that the experiments in [12] , where 
the exchange splitting could be seen, was conducted with 
a SWNT of L ~ 300 nm and a comparably small diame- 
ter of 0.8 nm. Similarly for samples A and B of pLlJ the 
reported lengths and diameters were La = 180 nm and 
Da = 1-1 nm and Lb = 500 nm and Db = 1-3 nm. 
In contrast, for the SWNT of sample C in [II], show- 
ing no measurable exchange effects, Lq — 750 nm and 
Dc = 2.7 nm were determined. For comparison, in the 
case of armchair tubes the diameters 0.8 nm and 2.7 nm 
would correspond to tubes with the wrapping indices (6, 6) 
and (20, 20), respectively. In the following we concentrate 
on large enough nanotubes and disregard any exchange 
effects. A detailed discussion of the processes leading to 
exchange splittings and their interplay with the bosonic 
excitations is postponed to a forthcoming article [15]. In 
the absence of exchange, a large degeneracy of the en- 
ergy spectrum is expected [16j, which in turn can be seen 
in a peculiar four-electron periodicity of the stability di- 
agrams (three equal in size Coulomb diamonds followed 
by a larger one) and of the Coulomb oscillation traces as 
discussed below. 

So far excitation lines of SWNT quantum dots were 
only addressed in [l3lfl4lfT7] in a meanfield approach. In 
the present work we do not only calculate the expected ex- 
citation lines but also give a quantitative calculation of the 
nonlinear current across a metallic SWNT quantum dot as 
a function of the gate and bias voltage beyond meanfield. 
As long as only the ground state energies are concerned, 
the meanfield treatment yields the right result but fails for 
a complete description of the excitation spectrum. Some 
results were already presented in abridged form in [Tg]. We 
focus on the low energy regime, which means that we con- 
sider the situation where only the gapless subbands of the 
SWNT are relevant. Since those bands exhibit to a good 
approximation a linear dispersion relation, the Tomonaga- 
Luttinger theory can be applied. For a typical SWNT this 



means a range of roughly 1 eV around the Fermi energy, 
which is still quite large if for example compared to the 
level spacing of 1.7 meV for a noninteracting SWNT of 1 
/im length. 

Starting from a microscopic description of noninteract- 
ing SWNTs with open boundary conditions we include the 
electron - electron interactions. The resulting Hamiltonian 
can be diagonalized using bosonization techniques |3J, ob- 
taining the spectrum and eigenstates of the isolated metal- 
lic finite size SWNT. We apply the so called constructive 
bosonization procedure jT8J (for a detailed review we re- 
fer to [H]), coping well with the discrete energy spectrum 
and especially with finite electron numbers. 

The electron dynamics of the quantum dot is obtained 
by solving the equation of motion for the reduced den- 
sity matrix (RDM) of the SWNT. The metal leads are 
described as equilibrated Fermi gases. Due to many-fold 
degenerate eigenstates of the SWNT Hamiltonian we find 
that offdiagonal elements (coherences) of the RDM can 
not be generally ignored. Apart from the degeneracy of 
the states, this is also a consequence of the interactions as 
we will show. The importance of taking into account co- 
herences for an interacting system weakly coupled to un- 
polarized leads was also recently discussed in [20], where 
transport through a metal grain is treated. The rates en- 
tering the master equation depend on the involved ener- 
gies and on the matrix elements of the electron operators 
in the SWNT eigenstate basis. Noteworthy is a variation 
of the tunnelling amplitudes along the nanotube axis, also 
if only ground states are involved in transport, depending 
on the complete energy spectrum of the collective excita- 
tions. 

The outline of this article is the following. In Section [2] 
we derive the master equation describing the dynamics of 
a generic quantum dot in second order of the tunnelling, 
gaining an expression for the current. In order to specify 
the explicit form of the master equation we derive the 
low energy Hamiltonian of metallic SWNTs in Section [3j 
and use the bosonization technique in order to diagonalize 
the SWNT Hamiltonian. The actual calculations for the 
current in the low and high bias regime are performed in 
Section |H 



2 Quantum dots 

In this section we show how the stationary current through 
the system described by the Hamiltonian JT]), see below, 
can be determined as a function of the electrochemical 
potentials in the leads and in the dot by calculating the 
dynamics of the RDM of the SWNT. Transport through 
the SWNT quantum dot can occur if the electrochemical 
potentials of source and drain are adjusted to different val- 
ues by the voltages V s and Vd- Since current can only flow 
if one transport direction is favored, the current-carrying 
system will have to be treated out of equilibrium beyond 
the linear response regime. Hence, the state of the quan- 
tum dot will in general be described by a nonequilibrium 
density matrix. For completeness and clarity, we show the 
determination of the corresponding equation of motion in 
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some detail mainly following [21]. The outcomes of this 
and of the following section will be used in Section |4] to 
obtain the I-Vb-V g characteristics of SWNT quantum dots. 

2.1 Model Hamiltonian 

In the following we examine the physics of a generic quan- 
tum dot, i.e. we consider a small metallic system weakly 
coupled to a source and a drain via tunnelling junctions. 
Moreover the electrochemical potential p g within the dot 
can be controlled by a capacitatively coupled gate voltage 
V g . We describe the overall system by the Hamiltonian 

H = H Q + H s + H d + H T + H g , (1) 

where Hq can describe an interacting SWNT (cf. equation 
l(44|) below) or another conductor with known many-body 
eigenstates. H s / d describe the isolated metallic source and 
drain contacts as a Fermi gas of noninteracting quasi- 
particles. Upon absorbing terms proportional to the ex- 
ternal source/ drain voltage V s / d , they read (I — s,d) 

oq 

where c qal creates a quasi-particle with spin a and energy 
£q,i = £q — eV s /d in lead s/d. The transfer of electrons 
between the leads and the central system is taken into 
account by the tunnelling Hamiltonian 

H T = E f d3r (Ti(rM(r)<P al (r)+h.c.) , (3) 

l=s,d a J 

where #£(r) an d ^L( r ) = E g ^( r ) C 9ffi are electron cre- 
ation operators in the dot and in lead /, respectively, and 
Ti{r) describes the transparency of the tunnelling contact 
at lead I. Finally, H g = —p g Af c = ~^cV g accounts for 
the gate voltage capacitively coupled to the dot, with M c 
counting the total electron number in the dot and c be- 
ing a conversion factor, that relates the electrochemical 
potential to the gate voltage. 

2.2 Dynamics of the reduced density matrix 

Our starting point is the Liouville equation for the time 
evolution of the density matrix p(t) of the system consist- 
ing of the leads and the dot. The tunnelling Hamiltonian 
Ht from equation (J3j) is treated as perturbation. We calcu- 
late the time dependence of p(t) in the interaction picture, 
i.e. we define 

p I (t) = U I (t,to) P (to)Uj(t,t ), (4) 
with the time evolution operator Uj(t,to), given by 

UAt t ) = e i{H & +H s +H d ){t-t ) e -j.{H & +H s +H d +H T )(t~t ) 

(5) 



with to being some reference time. Using (0} and §5§ the 
equation of motion becomes 

i»^=[^(t),/(t)], (6) 

Equivalently we can write 

P '(t) = /(to) - -J* [H^pih)] dt x . (7) 

Reinserting back into |(6|) yields 

pi{t) = - l -[H I T ,p{to)} 

+ (~) j*dh [H^t), [H^Hlp 1 ^)]] . (8) 

Since we are interested in the transport through the cen- 
tral system, it is sufficient to consider the RDM Pq of the 
dot, which can be obtained from p 1 by tracing out the 
lead degrees of freedom, i.e. 

Pq = Tr| eads //. (9) 

In general the leads can be considered as large systems 
compared to the dot. Besides we only consider the case 
of weak tunnelling, such that the influence of the cen- 
tral system on the leads is only marginal. Thus from now 
on we treat the leads as reservoirs which stay in thermal 
equilibrium and make the following ansatz to factorise the 
density matrix p 1 (t) of the total system as 

/O) = Pe(t)PsPd =■ /(pleads, (10) 

where p s and pd are time independent and given by the 
usual thermal equilibrium expression 

e -p{H 3/d -n B/d N' 3/d ) 

Ps/d = 7 , (11) 

^s/d 

with (3 the inverse temperature. As it can be formally 
shown [21], the factorization ((To]) corresponds, like Fermi's 
Golden Rule, to a second order treatment in the perturba- 
tion Ht- Furthermore, we can significantly simplify equa- 
tion (HJ) by making the so called Markov approximation 
(MA). The idea is that the dependence of p 1 it) on p 1 (f) 
is only local in time. In more detail, p 1 (t') is replaced by 
p T (t) in ([8]). The MA is closely connected to the correla- 
tion time r of the leads. In our case r is the time after 
which the correlation functions 

(^(r,T)^(^',0)) th := 

Tr le3ds ($l l (r,T)<P al {r',0)p le3ds ^ 

are vanishing and therefore erasing the memory of the 
system. If p 1 (t') is not considerably changing during r the 
MA is valid. It must be noted that the MA leads to an 
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averaging of the time evolution of on timescales of 

the order of r, such that details of the dynamics on short 
time scales are not accessible. Since we are interested in 
the dc current through the system, this imposes no restric- 
tion on our purpose. Finally we get, by plugging equations 
(J9j) , lfT0|) and fTTj) into (JSj) , the following expression for the 
equation of motion for the RDM, 

/>©(*) = -jT r \ sa ds [HT,Po(to) Pleads] 

+ (£\ Tr leads £ dti [H^(t), [f4(tl),p£,(t)pieads]] • 

(12) 

The first term vanishes because ( ( &al(i')) th — 0. Since we 
are only interested in the longterm behavior of the system 
we send to — > — oo. Writing out the double commutator in 
lfT2|) according to 

[A, [B, C}} = ABC + CBA - ACB - BCA, 

and introducing the variable t' = t — t\ one obtains 

i f°° 

- ^ Pleads / dt' [(H^(t)H^(t - <')P©(0Pleads + k.C.) 

- (H^p^PtedsH^t - t') + h.C.)] . (13) 

Now we insert the explicit form of Ht from equation ((3J) 
into lfl3|) and perform the trace over the lead degrees of 
freedom. Thereby we exploit the relations 




as well as 

=0 ioral^a'l', 

such that 

l=s,d a- J J J0 

{[£ al (x,y,t')^(x,t)V a {y,t - t')p G {t) 
+ T al (x,y,t')& a {x,t)^{y,t - t')p Q (t)] + h - c - 
- [T; i (x,y,t')^(x,t)p I Q (t)1r a (y,t-t') 
+ S^y^'^ix^p^t^Ky^-t'^-h.a}. (14) 

In lfl4|) we have introduced 
£ai(x,y,t') :=Ti(x)Tt(y)(^i{x)0U(y,-l/)) = 

\ I th 



with pp(e) = — where p/(e) is the density of 

energy levels in lead I and /(e) is the Fermi distribution. 
Similarly 

T al (x,y,t') :=T^x)Ti(y) ^{x^y.-t'))^ = 
T^xmy) Jdepf(e)Y,P q {*)<t> q {v)e Ue - eVl)t '- (16) 

Here pf{e) = pi{e)f(e). In order to proceed, it is conve- 
nient to represent the RDM in the eigenstate basis of the 
dot Hamiltonian Hq. Assuming that we can diagonalize 
the many-body Hamiltonian Hq (see Section [3]), this al- 
lows us to extract the t and t' dependence of the electron 
operators in lfl4|) . To proceed further, we carry out the 
following two major approximations: 

I. ) We assume that matrix elements between states 
representing different charge states vanish (the number 
of electrons in the dot influences the electrostatics of the 
whole circuit, hence is "measured" permanently). 

II. ) The so called secular approximation is applied, i.e. 
we only retain those terms of lf!4|) . which have no oscil- 
latory behavior in t. The latter causes that coherences 
between nondegenerate states are not present in the sta- 
tionary solution of (fT4|) . For the dynamics this means that 
we can not resolve the evolution of p Q (t) on time scales of 
h/(E m — E n ), where E m and E n are two distinct energy 
levels of Hq . 

In the end we can divide p@(t) into block matrices 
Pq En (t) , where En denotes an energy level of the isolated 
dot containing N electrons. To simplify the notation we 
give the resulting equations of motion in Bloch-Redfield 
form, 

kk' 

+ E EE^S^w. (17) 

M=iV±l E> kk' 

where the indices n, m, k, k! refer to the eigenstates of Hq. 
The Redfield tensors are given by (I = s, d) 

p£jv \ \ fir r ( + )E N E' M r r (-)E N E' M \ 

n-nmkk' — /-^ 2-~i \ Umk ' 1 l,njjk ' ° nk± l,k< jjm l> 

I M,E'.j 

(18) 

and 

nE N E' M _ V"^ r (a)E' M E N , 1Q v 
^nmkk' ~ 2-~i l,k'mnk ■ l la J 
l,a=± 

The quantities determine the transitions be- 

tween states with N particles and energy En to states 
with M particles and energy E' M . In detail we obtain for 
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transitions N — > N + 1 

(a)E N E' N + 1 _ _j_\T^ f ,3 
L Ik' rank ~ ^2 / 



Since the relation J", 



(+)Ejv -Ejv±i 



r, ( :]f wJ3Wi V holds, 



we arrive, by inserting (|24| into (|23| . at the following ex- 
pression, 



oM*))f"/" +i H(y));r iCN x 



dt' T a i{x,y,i!)e 



i' 



(20) 



ii = 

^2eRe 



N.E.E' 



r, 



(+)E N E' n 



l,njjk 



r, 



(+)Ejv E' N 



hnjjk 



) Pkr, 



(25) 



For transitions N — > A" — 1 is 

-,(a)E N Ejv-i 



/ k' rank 



mx)) E k ^(Mv))nr EN 



dt'£ rTl (x,y,t')e- a U E '»-i- E ») t ' , (21) 



where we have defined the matrix elements 

with the states \k) and |m) having energy En, E' n+1 and 
particle number N, N + 1 respectively. 

Equation ifTTl) governs the dynamics of the S WNT elec- 
trons. In the following we deduce therefrom the current 
through the system. 



For the actual calculations we are going to replace p ! < EN 
in (|25|) by the stationary solution p J s l N of (fT4| . because 
we are only interested in the longterm behavior of the 
system. 

Until now our treatment has been quite general, hav- 
ing made no assumptions about the nature of the dot so 
far. But the actual transport properties depend on the mi- 
croscopic structure of the system. As seen from equations 
(|20)l and l|2"T]l . we shall have to determine the spectrum 

(k \&l(x)\ m) , (22) of Hq and the matrix elements (^(x))^ En±1 i n order 

to solve master equation (fl7]l . This is the subject of the 
next section. Moreover the leads and the geometry of the 
tunnelling contacts will influence the system via the quan- 
tities T a i and E a i as discussed in Section l4~T1 and Appendix 

El 



2.3 Current 

The current is essentially the net tunnelling rate in a cer- 
tain direction at one of the leads. Thus the current at lead 
I will be of the form 



jv 



■N+l 



■N-l\ 



(23) 



where we use the convention I = s/d = ±1. On the long 
run the currents at the two leads have to be equal, oth- 
erwise charge would accumulate on the dot which is pre- 
vented by the charging energy. The rates jj^^ n±1 can 
be obtained from the time evolution of the occupation 

probabilities Pn — Tr (pqJ , where p Q ,N is the RDM 

for states containing A" electrons. In more detail, the oc- 
cupation probability of the charge state N is reduced by 
tunnelling events changing the number of electrons from 
A^ to N ± 1 and is increased by processes transferring the 
charge states N ± 1 to N, hence 

Pn = Tr fp^ N 

I 

Using Tr (pq N ) = Eb Tr (pq En ) together with (17]) we 
can easily identify the rates appearing in (f23|) : 



-iiV-*JV±l 



^(+)E N E N±1 I t E N . I ,E N - r (-) E N E N±1 



/A 1 l.niik Pkn "T" Pnk 1 l,kjjn 



E,E' nkj 



(24) 



3 Low energy description of metallic finite 
size SWNTs 

In the following we derive the Hamiltonian Hq describing 
the interacting electrons in a metallic finite size SWNT at 
low energies. Typically (depending on the diameter of the 
nanotube) "low energies" mean a range of 1 eV around 
the Fermi energy, such that the dispersion relation for 
the noninteracting electrons is linear [3] around the Fermi 
points. Since we are interested in the transport across a 
finite size system, open boundary conditions (OBCs) for 
the wave functions must be imposed and we can deduce 
the wave functions and the energy spectrum of the non- 
interacting system. Switching on interactions we still can 
diagonalize the Hamiltonian by using the so called con- 
structive bosonization procedure |18|ll9j , which we also 
employ in order to rewrite the electron operators in a suit- 
able form for the determination of the matrix elements 



3.1 Noninteracting SWNTs 

The bandstructure of SWNTs is conveniently derived from 
the one of the unbound 2p z electrons in graphene sheets 
[1]. Since each unit cell of the graphene lattice contains two 
atoms p — 1,2, the band structure consists of a valence 
and a conduction band which touch at the corner points 
of the first Brillouin zone. Two of those Fermi points, 
F = ±Kq, are independent, i.e. don't differ by a recipro- 
cal lattice vector. The energy dispersion around the Fermi 
points is very well linearly approximated with a Fermi ve- 
locity vf — 8.1 • 10 5 m/s. Since SWNTs are essentially 
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Fig. 2. On the right, energy spectrum of a SWNT with open 
boundary conditions (OBCs) described in terms of left (L) and 
right (7?) branches. It is constructed from suitable combina- 
tions of travelling waves (cf. equation {29}) whose spectrum is 
shown on the left side. 



graphene sheets rolled up to a cylinder of a certain diam- 
eter, we obtain the SWNT band structure by imposing 
periodic boundary conditions around the circumference 
L± of the tube, leading to the quantization of the allowed 
transverse wave vectors 

2tt 

k ± = —m, m = 0,±l,±2... , (26) 
L± 

and thus to the formation of subbands labelled by to. For 
metallic SWNTs the Fermi points F will satisfy condition 
f26|) , hence the corresponding valence and conduction sub- 
bands will have no gap. In the following the focus is on 
armchair SWNTs. Then only the gapless sub-bands with 
linear energy dispersion nearby the Fermi points are rel- 
evant [HE]- At each Fermi point there are two different 
branches r = R/L associated to right and left moving 
electrons. The corresponding Bloch waves are of the form 

<Pr/l,f,Jt) = e lK >fl/L,F(r-), (27) 

where K measures the distance from the Fermi points ±Kq 
(Fig. El left). 

In more detail, the Bloch waves ^Pr/l.f{ t ) a t the Fermi 
points are given by a superposition of wave functions living 
either on sublattice p — 1 or p — 2, 

Vr/lA*) = ~m= E elF ' R fprFX(r — R — t p ), (28) 

where Nl is the number of lattice sites which are iden- 
tified by the lattice vector R. The position of the atoms 
p = 1, 2 in the unit cell is given by t p and the func- 
tions x are the p z orbitals. The value of the coefficients 
fprF depends on the considered SWNT type. For sim- 
plicity we concentrate on armchair SWNTs for which we 
have fi r F = 1/V2, farF = — 1/ v / 2sgn(rF) , where we use 
the convention that R/L = ±1. In this article, we are 
interested in finite size effects and therefore we impose 
OBCs instead of periodic boundary conditions on the sin- 
gle electron wave functions. Generalizing [22] to the case 
of SWNTs we introduce standing waves which fulfil the 
OBCs (Fig. H right): 

^ k/L,S r ) = [<PR/L,K ,*( r ) ~ ( PL/R,-K ,- K (r)] , 

(29) 



with quantization condition 

k = 7t(to k + A)/ L, m K = 0,±1,±2,..., 

where L is the SWNT length. The offset parameter A 
occurs if there is no integer n with Kq = (nn/L)e\\ (where 
e II is the unit vector along the tube axis) and is responsible 

for a possible energy mismatch between the states of the R 
and L branches, defined by relation (|29| . From the linear 
energy dispersion relation around the Fermi points, the 
energies of the standing waves follow as 

s fK = sga(f)Hv F K. (30) 

Including the spin degree of freedom, the electron op- 
erator reads 

f=R,l K ' a a 

where the operator Cf aK annihilates the state \f^ c ) \<r). 
Since the Bloch waves can be divided into a slowly and 
a fast oscillating part (cf. ([27)1 ) we can as well split off 
a slowly varying part from the operators ^(r), i.e. we 
introduce the ID operators 

^„f{x) = E e lsgn(F)Ka: c, ffK , (32) 
in terms of which the 3D electron operators read 

<^(r) = %/I^sgn(F)x 

F 

[<Psgn(F)R,F(r)lp Rr7F (x) + lfis g n(F)L,F(r)lpLaF( X )] ■ ( 33 ) 

Later on, the bosonizability of t/wf^) be of quite 
some significance for the calculation of the transport prop- 
erties. 

From f30|) the appropriate Hamiltonian for the nonin- 
teracting system is easily derived. It reads 

H Qy0 = hv F sgn(r) ^ nc\ aK c faK . (34) 

fa k, 

3.2 Interacting SWNTs 

For the inclusion of the electron - electron interactions we 
have to add the following term to the SWNT Hamiltonian, 

\Hl d * r j d 3 r'^ a (r)^l,(r')V(r-r')^(r')K(r), 

(35) 

where V(r — r') is the possibly screened Coulomb poten- 
tial. Upon inserting lj33j) into lf35j) . integration over the 
coordinates perpendicular to the tube axis yields the in- 
teracting Hamiltonian expressed in terms of ID operators 
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Fig. 3. Scheme of forward scattering terms. All those processes 
conserve the number of electrons in each branch. 



and an effective ID interaction V e s(x, x'), see (|37| below. 
As in p] we expand the Bloch waves into their sublattice 
contributions (|28|) and ignore the difference between the 
intra and inter lattice correlations, which are only of rel- 
evance at the length scale of the next neighbor spacing 
of the carbon atoms. The final form of Vq we obtain by 
ignoring all fast oscillating terms in the interaction, or in 
other words by keeping only the so called forward scatter- 
ing processes associated to the effective potential V e g as 
depicted in Fig. [31 This leads to an expression of Vq in 
terms of the ID densities pf a F{x) — ^l aF (x)^ 7 ~crF (x) , 



V, 



Jeee 



FF' 



dxdx' pf vF (x) V e ff (X , x')pf 'v'F' («' ) • 



(36) 



The approximations that have been made for getting from 
|3o|) to ([Ml mainly mean that we neglect, as already men- 
tioned in the introduction, any kind of exchange effects 
and thus are only valid for large enough SWNTs. Using 
|28|) and l[3"Tj) the effective potential V c h is obtained from 
an integral over the coordinates perpendicular to the tube 
axis of the 3D Coulomb interaction weighted by the p z 
orbitals, i.e. 

V e g(x,x') = ^x 

E / d 2 r ± f d 2 r' ± \ x (r - R)\ 2 F(r-r') | X (r' - R')f . 

R,R' 

(37) 



3.3 Bosonization 

Here we show how the introduction of bosonic excita- 
tions enables us to diagonalize the SWNT Hamiltonian 
Hq = Hq$ + Vq by recasting it into a sum of a fermionic 
and a bosonic part. Besides we give the identity that ex- 
presses the electron operators in terms of the boson op- 
erators. The connection between fermionic and bosonic 
operators can be obtained from the Fourier coefficients of 
the electron density operators. To be more specific, we 
Fourier expand the electron density operator, 



PfaF(x) 



— V 

2L ^ 



„isgn(F)qx 



(38) 



with q = fjn g , n q G Z and define the operators 
1 



^crsgn(f)g 



zPfasgn(r)qi Q ^> 0. 



(39) 



The operators fulfil the canonical bosonic commuta- 



tion relations 



baq,bl, q , 



5aa'5 qq i as it can be shown 



[19| by using the explicit expression 

basgn(r)q — __ c f,j K c faK+sgn(f)qi q > 0- (40) 



From the previous equation we can see that the 6, 6^ op- 
erators annihilate or create collective particle hole excita- 
tions within the branch {fa). 



3.3.1 Diagonalization of the SWNT Hamiltonian 

As for example shown in [19] for a generic ID system, the 
noninteracting part of the Hamiltonian is already diagonal 
in the bosonic operators. For the SWNTs in particular, we 
obtain from (021) and lt40l). 



0,0 



so 



E 



Ekl^A? + E 

q^O f 



-M?„ + Asgn(f)Af fl 



(41) 



where £o = fopf is the level spacing of the noninteracting 
nanotube. The first term in (fUJl describes collective par- 
ticle hole excitations, whereas the second term represents 
the energy cost of the shell filling due to Pauli's princi- 
ple. Specifically, A/f CT = Y1 K c lcrn c f<?n counts the number 
of electrons Nf a in the (fer)-branch. Without loss of gen- 
erality we have assumed that we only have states with a 
positive number of particles in each branch. 

Plugging the Fourier expansion of the density opera- 
tors ([381 into ([36"1) and taking into account the definition 
of the 6-operators yields the bosonized form of Vq , 



2 •-</■'??' 
<?>o 



EE (^sgn(f)g + ^Lgn(f)<}) (^o-'sgn(f ')q + ^t'sgn(f')g) 



<7<7' rr 1 



(42) 

Here E c — Wqq is the SWNT charging energy responsi- 
ble for the Coulomb blockade and the effective interaction 
potential is absorbed into 

W qq = -^- j dx ( dx'Vcfi (x, x') cos(qx) cos(qx'). (43) 
■Li Jq Jo 

Since the bosonic operators appear only quadratically in 
(|4l"l) and (J42j) , the Bogoliubov transformation [23] can be 
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applied to diagonalize the bosonic part of the total SWNT 
Hamiltonian. Here we only give the result, namely 



1 



EM 2 



so 



TO 



J r<7 

2 



Asgn(f)Aff c 



q>0 j=c,s S=± 



(44) 



The first line of l144[) describes the energy cost to add new 
particles to the system. The excitations are created by 
the bosonic operators o}- Sq . Four channels are associated 
to total (j5 = c+, s+) and relative (jS — c— , s— ) (with 
respect to the occupation of the R and L branch) charge 
and spin excitations. To be more precise, the new opera- 
tors djSq are related to the old operators & CTS gn(r) g via the 
Bogoliubov transformation: 



J f7sgn(f)g 



(45) 

jS 




jS = c+, c— , s+, s 
fa = RtR 1,11,1 4' 

(46) 

For the transformation indices Sjs q and Cjg q we get in the 
case of the three modes jd = c— , s+, s— , 

Sjsq = 1 and C jSq = 0. (47) 

c+ there is an interaction dependence, 




Generalized spin-charge separation occurs, since for the 
three interaction independent channels j5 = c— ,s+,s— 
the energy dispersion is the same as for the noninteracting 
system, 



Ej&q 



HvfQ = hvFyriq 



'Oq, 



1,2, 



but the energies of the c+ channel are enhanced by the 
repulsive interaction, 



£c+q 



e q(l + 8Wqq/e ) 1/2 . 



(49) 



The ratio between e c + q and So q for a (20,20)-armchair 
SWNT is shown in Fig. (j4]). For the actual calculation 
of e c +q we have substituted |x(r — R)\ 2 by S(r — R) in 
(|37j) . The extension of the p z orbitals has been modeled 
by introducing the average distance a± ps 0.15 nm of the 
electrons from their nuclei [31, i-e. we used for the 3D in- 
teraction potential the following expression 



V(R-R') 



\R-R' \/L, cr , 



47re e 



R-R> 




300 



Fig. 4. The ratio e c +q/so q as a function of Eo q for a (20,20) 
armchair SWNT of 980 nm length. Here we show the decay of 
£c+ q/so q for an unscreened and a screened (screening length 25 
nm) Coulomb interaction. In both dielectric constant 

e of 1.4, see Ref. [3], is assumed. 

Here ^screen is the screening length of the potential and 
eos is the dielectric constant. Due to the finite range of 
the interaction potential we find a decay of e c + q /so q with 
increasing q, which should be taken into account if higher 
excitations in the c+ channel are involved in transport. In 
the Luttinger liquid theory e c + q /eo q is usually assumed 
to be equal to the constant 1/g, where 



g = (1 + 8W qq /e y 1/2 , 



q = T 



(50) 



The eigenstates of the SWNT Hamiltonian in lj44"ll are 
\N,m):= [] (rrWr 1/2 (a] Sq ) mjSq \N,0), (51) 

q>0,jS 

where | N, 0) has no bosonic excitations and the vector N 
defines the number of electrons in each of the four branches 
(fa). 



3.3.2 Bosonization of the electron operators 

The determination of the transport properties through the 
SWNT quantum dot involves the calculation of the matrix 
elements of the electron operators between the eigenstates 
JUJ. Therefor we give the so called bosonization identity 
for the operators fepfi), 



lpraF{x) 



1 



(52) 

The way (|52| is derived can be found e.g. in [19]. Here 
a is a convergence factor needed to ensure that only ex- 
citations with physically sensible energies are considered. 
It will be set to zero at the end of the calculation. The 
operator rj fa is the so called Klein factor. Its main effect 
if acting on a state \N,m) is to annihilate a particle in 
the (fer)-branch, more precisely 



\N.m) 



I AT- 
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where we use the convention j = R [,L \,L J, = 
1,2,3,4. The fact that the Klein factor annihilates fermions 

after all, expresses itself in the factor (— l)^--j= 1 Nj . Kf aF (x) 
yields a phase depending on the number of electrons in the 
(fa) branch, in our case, 



K faF (x) = 



J^sgn(F)(sgn(f)J\r ftr +A)x 



2L 



At last we have the fields (j) raF (x) given in terms of the 
bosonic operators b aq , 



i(j> raF (x) 



E 

q>0 



,-aq/2 



isgn(fF)qx h 



'(Tsgn(f)q- 



3.3.3 The matrix elements of the electron operators 

In the following we are going to determine the matrix el- 
ements (& a (x))f£, En±1 ) appearing in the transition rates 



1(20]) and (J2U) . As we know from equation ((33)) the 3D elec- 
tron operators "^(cc) can be written in terms of the ID 
operators Tp r aF{x) yielding 



i^)) E k k E ' N±1 =VI£sgn(F)x 

F 

tP Sgn (F)R,F( r ) (^RaF^fkk' ^ 
+ Vs g n(F)L, F (r) {tpl aF {x)) E ™, 



(53) 



Now the bosonization procedure used for the diagonal- 
ization of H & and the redefinition of ip r <rF(x) in terms of 
boson operators pays off. Denoting the SWNT eigenstates 
\k) and \k') as \N,m) and \N',m') according to equation 
(joTj) . we get, as shown in Appendix [Bj 



(N,m\ip fr7F (x) \N',m') = S N+e ,^ N ,(-l)^7=i 1 

e _ ^5I< J >o e aq J2js\ x f S J , F( x )\ 
Kf F (N'),A x ) 7 x 



N, 



= :A{x) 

IHI^WvfW^A.^)' (54) 

«>0 jS 

The parameters \^ F (x) (cf. Appendix [C]) for the three 
neutral modes jS = c— , s+, s— are given by 



H aF ( 3 



^isgn{Fr)qx 



-A J - 



(55) 



1 

0.8 
0.6 
0.4 





L = oo 




~ ~ " screen 




L = 25 nm " 




screen 








> 


f * 

















0.2 



0.4 0.6 

x/L 



o.s 



Fig. 5. The non-oscillatory position dependence of the matrix 
elements {k \ipf-crF{x)\ k') is described by A(x), if no c+ excita- 
tions are involved. For a repelling interaction, A(x) is strongly 
suppressed at the SWNT ends. Here we used the values of 
Scq/eoq of a (20,20) armchair SWNT (L = 980 nm) as shown 
in Fig. (J4|) . 



The function F(X,m,m'), describing the dependence of 
the matrix elements on the bosonic excitations, can be 
expressed in terms of the Laguerre polynomials L™ as 
we show explicitly in Appendix [Bj see equations J93|) and 
([94)) there. It is interesting to note that the interaction 
leads to the formation of a non-oscillatory x dependence of 
the matrix elements (k \ipfaF(x) \ k') . For matrix elements 
between states with no bosonic c+ excitations this effect 
is described solely by the function A(x). From expressions 
(|55| and (|56| for the parameters \^ F we find 



A(x) 



-(f^cos^ + f^sin 2 ^)-!) 



(57) 



and for the c+ mode we have 



Since we are considering a finite ranged interaction, the 
ratio £ C q/soq goes to 1 for large q and so the sum in ([57I 
converges even for a = (using a finite a in order to dis- 
regard unphysical states outside of the low energy regime 
has only little effect on A(x)). In Fig. [5] we show A(x) for 
a (20,20) armchair SWNT with screened and unscreened 
Coulomb potential, using the values of £ C q/£oq as shown in 
Fig. Q. It is evident that A(x) (and hence the tunnelling 
amplitude at low energies) is smallest at the tube ends and 
increases towards the middle of the nanotube. For an at- 
tractive interaction the opposite behavior would be found 
and for the noninteracting system A(x) = 1. Hence A(x) is 
closely related to the dispersion relation for the excitation 
energies e cq . However, in this article the contact geometry 
is fixed and so the actual form of A(x) doesn't influence 
the transport properties qualitatively, as we will show be- 
low when discussing the role of the tunnelling contacts. 



4 Current through a SWNT quantum dot 



\f£ F (x) — 

1 



' cos(qx) + i , 

£0q 



' -^-sgn(Ff) sin(gx) 

£ 



Knowing the spectrum of H Q and the expressions for the 
electron operator matrix elements we can return to Sec- 
(56) tion 12.21 and perform the actual calculation of the current 
through a SWNT quantum dot. 
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4.1 Influence of the tunnelling contacts 

Due to the integration J d 3 x J d 3 y, the rates from (TSUI and 
lf2"Tj) depend on the geometry of the tunnelling contact. At 
low enough energies (i.e. as long as no degenerate states 
with very different bosonic excitations or fermionic con- 
figurations are considered), the matrix elements of the ID 
operators are slowly varying compared to the extension 
of the tunnelling contacts described by the transparencies 
Ti(r). Then the influence of the tunnelling contacts at the 
tube ends can be taken into account by introducing the 
parameters (cf. Appendix |A]) 

$iff>(e) = 

J d 3 r J d\'Tt{r)T l (r')Y,^ q {r)4> lq {r')A{x)A{x') 
x sgn(FF')(p sgn{F)r ^ F (r)(p* sgn{F , )r , F ,(r')iyi FF ,(A), 

FF' 

(58) 

where the phase factor viff'{A) = e lS '- d ™ sn( - F - F ' )A ac- 
counts for the band mismatch A. With the new parame- 
ters <Pjff(e), equation pfjj) can be rewritten as 

(a)E N E' N + l _ 
I k'mnk 

pE/ dspf(e)<P lff ,(e) (A*lf k r m E ' N+1 X 

rr'a 

(4'S' N+lEN r ti>e a H*-' v *-&™- E »)y (59) 

V r a J nk J a 

with pf(e) = pi(e)f(e) as defined below equation (fl6|) . 
Additionally we have introduced the notations ipfui := 
ipr<rK (xi)/ A(xi) and xi — 0, L for I — s, d. For l[2l"j) the 
parameterization can also be performed yielding, 

p(a)E N Ejv-i _ 
I k'mnk 

rr'a 

(Mnk- lEN / cft'e^H+^-K-!-^))*', (60) 
Jo 

where pf{e) = pi{e){l — /(e)). The integrals over e and t' 
in (|59ll and (|60|) can be carried out by using 

J deg{e) J dt' e a ^ £ - E)t ' = nhg(E)+aihV J j^de, 

(61) 

where V denotes the Cauchy principal value. The first 
term on the right hand side of 1|61|) corresponds to pro- 
cesses which conserve the energy. Additionally, we have 
the principal value terms which can be attributed to so 
called virtual transitions since they cancel in the expres- 
sion for the current but nevertheless can affect the trans- 
port properties indirectly via the time evolution of the 
RDM. 



From now on we assume furthermore that the leads 
are properly described by a 3D electron gas (think of gold 
leads for example) . In Appendix [A] we then find for a re- 
alistic range of the lead electron wavenumbers q, 

$lr?'{e) = (e), (62) 

with 

^-^iTKfi + Tp)! 2 ^ 2 ^), (63) 
a,p 

simplifying equations lf59"ll and l|6"0l) further. We should 
note that the case 

$lf?>(e) 8fr$i(e) (64) 

corresponds to leads which are polarized with respect to 
the band degree of freedom f , in complete analogy to spin 
polarized leads. Hence it will be interesting to discuss in 
the future for which type of contacts, (|64| eventually could 
be fulfilled. 



4.2 Excitation lines 

If certain transitions are possible depends among other 
things on the available energy and hence on the applied 
gate and bias voltages. From lj59j) and ([60|) together with 
the evaluation of the appearing L dt' . . . integral accord- 
ing to ([6lj) , we know the resonance condition for tunnelling 
in/out of lead I, In detail, at temperature T = 0, transi- 
tions from a state with N electrons and eigenenergy En 
to a state with M = N ± 1 electrons and eigenenergy Em 
are possible under the condition 

eVi < E N - E N+1 + fig, for M = N + 1, 

eVi > E N -i - E N + fig, for M = N - 1. (65) 

We expect that for low temperatures the current can only 
change considerably at the position of the correspond- 
ing lines in the eV - fi g plane. But we should mention 
that in principle the so called virtual transitions (cf. (|6T|l ) 
are sensitive to the values of the bias and gate voltage 
outside of the excitation lines, if coherences of the RDM 
influence transport. However, in our case of unpolarized 
leads we do not find a significant change of the current 
between two excitation lines. Furthermore not all of the 
resonance conditions lead to a considerable change of the 
current at the corresponding lines in the eVf, - fi g plane. 
Relevant transitions are those between states with consid- 
erable overlap matrix elements {N ,m\ij:f aF {x)\N' ,m') 
and for which the occupation probability of the initial 
state is large enough. 

In order to add another electron to a N particle ground 
state the extra energy 8fiM — E N+1 — E N — (E N — E N _ 1 ) 
has to be payed. Here E M denotes the ground state en- 
ergies. From the SWNT Hamiltonian H@, equation ([44]) . 
we deduce 

Sfi4m+2 = E c + 2\A\e , 

6fi 4m = E c + (l-2\A\)s . (66) 
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The energy S[1n is a direct measure of the height and the 
width of the Coulomb diamonds in the eVb - p, a plane. 
Thus a repeated pattern of one large Coulomb diamond 
followed by three smaller ones is expected for A = 0, 1/2. 
Otherwise the pattern will consist of a large diamond fol- 
lowed by a small, a medium and again a small one. In the 
experiments of [TT] sample C showed the first pattern re- 
peating very regularly, whereas samples A and B revealed 
the second pattern. 

For sample C of |1 lj not only the Coulomb diamonds 
but also a bunch of excitation lines could be resolved. In 
[TT] the positions of the Coulomb diamonds and of the low- 
est lying excitation lines were determined. The mean field 
theory of [13] was used for comparison. Apart from the 
height of the large Coulomb diamonds, all the lines from 
the experimental I — V characteristics of sample C could 
be reproduced by an appropriate choice of five mean field 
parameters E c , eo, A, J and dU. The first three parame- 
ters also appear in our theory. The mean field parameters 
J and dU are the exchange energies with respect to the 
spin and band degree of freedom. For sample C the choice 
of pi] was E c = 6.6 meV, e Q = 8.7 meV, 2Ae Q = J = 2.9 
meV and dU — meV. With our theory exactly the same 
Coulomb diamonds and excitation lines are recovered by 
choosing E c — 9.5 meV, £o = 2.9 meV and A = 0. We 
think that in this case our choice of parameters is much 
more realistic than the one made in [TT], with an unreason- 
ably high J of 2.9 meV (compare this to J = 0.5 meV in 
[12J, where a considerably smaller SWNT was examined). 
We therefore conclude that our treatment of the interac- 
tion, where only forward scattering events are considered 
and exchange contributions are not present, is valid here. 
In Section [4~4l Figure [9] shows the numerical result for the 
corresponding current. 



4.3 Low Bias Regime 



N=4m 



A=0: 



L R 

* L * L 

i * i if 

HH 

(lx) 



A=l/8: i_ 



(lx) 



N=4m+1 


N=4m+2 


N=4m+3 


L R 


L R 


L R 
▲ i ▲ 


i_ L i__L 
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tin 


i_ L i__L 
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(4x) 


(6x) 


(4x) 












tttt 


(2x) 


(lx) 


(2x) 



Fig. 6. Ground states of the SWNT for different numbers of 
electrons. For each charge state we show one representative 
ground state, on top for aligned bands and on the bottom 
for mismatched bands. In brackets the corresponding degener- 
acy is given. Full black and orange arrows represent occupied 
states. The bold (orange) arrows indicate electrons which can 
contribute to transitions N — > N — 1 in the low bias regime. 
Incoming electrons for the transitions N — > N + 1 can be ac- 
commodated by the states represented by the dashed arrows. 
So the number of bold and dashed arrows is equal to Cn,n-i 
and Cn,n+i respectively. 



of the groundstates containing M particles will all be the 
same in the stationary solution. Therefore, we introduce 
the probability for finding the system in charge state M 

by 

Pm = d M pNN{t), 

where g?m is the degeneracy of the corresponding ground 
states. Using lfT7|). the general expression for the master 
equation, we find the following CME for Pm, 



In the following we consider the low bias regime, i.e. the 
bias voltages and the temperature are low enough that 
only ground states with N and N + 1 particles and ener- 
gies Em, Epf +1 can have a considerable occupation prob- 
ability. In this case the importance of taking into account 
the off-diagonal elements of the RDM in (fTF)) depends cru- 
cially on the parameters 4>i rr '(e) from l(58|) . From the ex- 
pressions l|59p and f60|) for the rates it is evident that for 
our assumption of unpolarized leads, i.e. under condition 
(|62|). the time evolution of the RDM elements between 
states with the same band filling vector TV is decoupled 
from elements between states with different TV. Since the 
current only depends on the time derivative of the diag- 
onal elements of the RDM, the elements mixing states 
with different N will have no influence on the current and 
therefore can be ignored. Because all considered ground 
states do have a different N, in the low bias regime we 
only have to take into account diagonal matrix elements 
in the master equation. We will refer to this kind of master 
equations as "commonly used master equations" (CMEs). 

In the low bias regime the occupation probabilities P^n 



P 



M 



-R r 



, ~ L NNNN j 

(1m um 



Pm r 



T?0 

M a M' 

NNN'N 



P 



M' 



N' 



til/' 



(67) 



With equations (fT8| and (fT9|) for the Redfield tensors we 
obtain 



P. 



M 



J l ~ 



d M r E° M ,E° M 

dM> 



2 EE ( r mw' p M - 5^ r ^» M PM' ) , (68) 

I N' ^ 



where we have defined 



\- p E° M E° 
N' 



r"M" M ' r (.+) E M E ° M i 

• i m -w- ■— 1 INN'N'N ■ 
N> 



Using the relation 



l JL sr r E °M; E M = \- r^-fu 



(69) 



(70) 



N' 



N ' 
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the stationary solution of lf68| is easily obtained, 



p _ t-~dN' I NN' 



with 



r to t - y^y^ ( r t £ N ? +1 + r t 



l N' 



7° J?0 7?0 
I NN' 



(71) 



(72) 



The current is obtained from l(23j) . Evaluated at the source 
for example it yields, 



In. 



N+l 



^N^N+l p^N + l^N 

N' I * s NN' d NN' 



P P TTlO pO 

1 s NN' 1 d NN' 



lot 



The rates fj are obtained by using equations ((59 

and (60]), 



1 I NN' 



>/( £ oe(^)^ +i (^);:; 



and 



r-IJ pu 



r 

z l NN 



NN' 



i -Q J4 1\ 



(,/> .^w^w + i 



Here we have assumed that and p/ are constant in the 
relevant energy range. Furthermore we have defined e; = 
eVi - AE and AE = E% - E° N+1 . From (54]) we obtain 
(remember that ^f„i := i>r<xK { x l) I A(xi)), 



fa N' 



where Cm,m' m the number of ground states with M' par- 
ticles that differ from a given band filling vector N for one 
of the ground state with M particles only by a unit vector 
(see Fig. [6]). Therefore we get 



T7.U pU 

I NN' = 7lj{£l)(;N,N+l 



(74) 



N' 



as well as 



2 E = 7/ (1 - /(ei)) Civ+i.iv, (75) 



N' 



with 7; = j^>ipi. Inserting the rates (74]) and ([75| into 
expression (75]) for the current results in 



l N.N+l — e 



Cn,n+iCn+i,n1s1cI [/(g s ) — /(grf)] 
7* [/(ei)CW,JV+i + (1 - Hei))C N+1 , N ] ' 



(76) 



4.3.1 Linear conductance 

In the regime \eV\ <C kT <C £o we can further simplify 
(|76| by linearising Inn+i m the bias voltage. We choose 
-eV s = eV d - ^ 



e^f and obtain 



lN,N+l 



e 2 f3CN,N+\CN+l,Nlsld 



-AE 



:V b , 



Ei H {foC N ,N+i + (1 - fo)C N+hN ) [ e ~ AE + 1) : 

with fo — f(—AE). Unlike one could expect, the maxima 
of the conductance Gn,n+i = lN,N+i/Vb are not at AE = 
0, but at 



AE„ 



1 Cn+i.n 
-ln- 



2(3 Cn,n+i ' 

which is only zero for Cn+i.n = Cn,n+i- The height of 
(73) the conductance peaks is 



-/V../V+1 



7s7dCjV,JV+lCiV+l,iV 



(7s + 7d) (CiV,iV+l + CjV+l,JV + 2y/CN,N+lCN+l,N) 



e 2 /?. 



(77) 



We still have to determine the values of Cn,n+i and Cn+i.n, 
which depend on the mismatch between the R and the L 
band. If both bands are aligned (A — 0, 1/2) one finds 
from Fig. H CW,iv+i = 4,3,2,1 and C n +i,n = 1,2,3,4 
for N — Am, Am + 1, Am + 2, Am + 3. Then the conduc- 
tance Gn,n+i shows fourfold electron periodicity with two 
equally high central peaks for N — Am + 1, Am + 2 and 
two smaller ones for N = Am, Am + 3 (cf. Fig. [7] a)). 
The relative height between central and outer peaks is 
GZi+i Am +2/GZ?A m +i - 27/(10 + 4V6) « 1.36. Note 
that this ratio is independent of a possible asymmetry 
7s 7^ "id in the lead contacts. In addition, the conductance 
is symmetric under an exchange of the sign of the bias 
voltage. 

In the case of an energy mismatch between the R and 
the L band exceeding well the thermal energy, the de- 
generacy of the ground states is either 1 or 2 and we 
get CVat +1 = 2,1,2,1 and C N ,N+i = 1,2,1,2 for N = 
Am, Am + 1, Am + 2, Am + 3. Therefore, according to (77]). 
all the conductance peaks have the same height as de- 
picted in Fig. [8] a) . 



4.3.2 Low bias regime well outside of the Coulomb 
diamonds 

Now we examine the regime where still only ground states 
are occupied but where we are well outside the region 
of Coulomb blockade, i.e. e > \eVi ± AE\ > kT. Then 
we have |/(e s ) - f(e d )\ = 1. If e.g. eV s - AE < and 
eVd — AE > 0, such that electrons tunnel in from the 
source and tunnel out at the drain, we have f(e s ) = 1 and 
f( e d) — such that (|76| becomes 



lN,N+l 



Cjv,jv+iCiv+i,Jv7s7d 

1sCn,N+1 + "fdCN+l.N 



(78) 
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Fig. 7. Gate traces for A = 0, i.e. for aligned R and L bands, 
at low bias voltage. Parameters are E c — 9.5 meV, ksT — 
0.10 meV, £q — 2.9 meV and the chosen asymmetry is 7 S = 
5jd = 4.9 • 10 10 s _1 . a) Conductance in the linear regime eVt <C 
ksT <C eo. Despite asymmetric contacts, we find the repeating 
pattern of two small outer peaks and two large central peaks, b) 
Current in the regime fcsT <C |eV5| <C £o- Asymmetry effects 
appear. The fourfold periodicity is retained. The upper and 
lower pattern correspond to opposite values of the bias voltage. 



The height of the plateaus in the current of Figures [7] b) 
and [U b) are described by (|78|) . If the R band is aligned 
with the L band, we still find a fourfold electron period- 
icity. But only for 7 S = jd the pattern with two central 
peaks and two smaller outer peaks is preserved. The cor- 
responding ratio of the heights is 3/2. If 7 S jd this latter 
symmetry is lost. 

For mismatched bands and 7 S = jd we find like for 
the conductance peaks that all current maxima are of the 
same size. In the case of asymmetric tunnelling contacts, 
7 S 7^ jd, a pattern of alternating small and large peaks is 
found. 

If we invert the sign of the bias voltage, the current is 
obtained by flipping its direction and exchanging j s with 
7d in (|78|) . Then, if 7 S ^ jd, the current does not only 
change its sign but also changes its magnitude because of 
Cn,n+i 7^ Cn+i,n- 



4.4 High bias regime 

In the bias regime eV > eo not only the ground states will 
contribute to transport but also states with bosonic excita- 
tions and band filling configurations TV different from the 
ground state configurations. We refer to the latter type 
of excitations as fermionic excitations. Since the number 
of relevant states increases rapidly with increasing bias 
voltage, an analytical treatment is not possible any more 
and we have to resort to numerical methods in order to 
calculate the stationary solution of the master equation 
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Fig. 8. Gate traces for A = 0.2, i.e. for mismatched R and 
L bands, at low bias voltage. Other parameters are E c = 9.5 
meV, ksT — 0.10 meV, eo = 2.9 meV and the chosen asymme- 
try is 7s = 5jd = 4.9 • 10 10 s _1 . a) Conductance in the linear 
regime eVj, <C k B T <C e - The conductance peaks are all of 
the same height, b) Current in the regime ksT <g \eVi\ <ti eo- 
The contact asymmetry leads to alternating large and small 
current maxima. 



(fl7| and the respective current. From (|65| we know that 
at low temperatures the current only changes consider- 
ably near the excitation lines given therein. Therefore we 
can reduce drastically the number of [eV, ^ g ) points for 
which we actually perform the numerical calculations, sav- 
ing computing time. In Figs.[9]a) and[10l the current as a 
function of the applied bias voltage and the electrochemi- 
cal potential in the dot is depicted. The chosen parameters 
for E c , eo and A are the ones we have obtained for fitting 
the data of sample C and sample A of [H], respectively. 
Hence Figs. l9l and fTOl show the current for SWNTs with 
the R and L band being aligned (A — 0) and mismatched 
(A « 0.17). In both cases a symmetric coupling to the 
leads, i.e. 7 S = Jd, is assumed. In the transport calcula- 
tions the lowest lying c+ excitations are taken into account 
using g = 0.21 for the Luttinger parameter (cf. equation 
([50)1 for the definition of g) . 



In addition we have also determined the current us- 
ing the CME, hence ignoring any coherences in the RDM. 
For the current corresponding to Fig. [9l the quantitative 
difference between the calculations with and without co- 
herences is considerable in the region of intermediate bias 
voltage as we show in Fig. [9] b) . On the other hand the 
deviation of the CME result from the calculation includ- 
ing coherences is by far less pronounced for the parameter 
choice of Fig. [TUl The crucial point here is that the charg- 
ing energy E c is smaller than eo, the level spacing of the 
neutral system. Then the subsequent considerations are 
not strictly valid. Now we explain why coherences can't 
be generally ignored if considering interacting electrons in 
aSWNT. 
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Fig. 9. a) Current in a bias voltage - electrochemical potential 
plane for the symmetric contacts case, b) Difference plot of the 
current with and without coherences. Here ksT = 0.01 meV 
and the interaction parameter is g ~ 0.21. Other parameters 
are as for Fig. [7] 
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Fig. 11. a) Scheme of relevant energy levels, if transitions from 
the ground state with N = 4m particles to the first excited 
state with N — 4m + 1 particles are energetically allowed, 
while no transitions from the iV = 4m + 1 ground states to 
excited states with N = 4m electrons are possible, b) Possible 
transitions between the N = 4m and iV = 4m + 1 electron 
states of the noninteracting system, which are energetically 
allowed in situation a). The degeneracy of the eigenstates is 
given in brackets. 
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Fig. 10. Current as a function of the bias voltage and the 
electrochemical potential in the SWNT for symmetric con- 
tacts. The parameters here are chosen to fit the positions of 
the Coulomb diamonds of sample A in [11], i.e. E c = 4.7 meV, 
2Aeo = 2.8 meV, £o = 8.2 meV. 

4.4.1 Why and when are coherences needed? 

As in the low bias regime we assume to have unpolarized 
leads, i.e. we use condition lj62j) . such that we can ignore 
coherences between states with different fermionic con- 
figurations N. Unlike in the low bias regime, we are not 
only left with diagonal elements of the RDM but there still 
might be coherences between degenerate states which have 
the same N but different bosonic excitations m. For the 
importance of these kind of coherences it is illuminating 
to discuss our system without electron-electron interac- 
tions, i.e. for the moment let us assume that an eigenbasis 
of Hq is given by the Slater determinants of the single 
electron states ) \a) . Furthermore, we concentrate 

without loss of generality on the case A = 0, and we as- 
sume that the charging energy E c exceeds Eq. Each of the 
Slater determinants can be denoted by the occupation n 
of the single electron states. In the case condition lf62]) 
holds, it is again easy to show that coherences vanish in 
the stationary solution of the master equation (fl7|) , if the 
RDM is expressed in the \n) basis. But of course we still 
could use the states |7V, m) from i|5ip as eigenbasis, now 
with four neutral modes c+,c—,s+,s— . In the \N,m) 



basis it is crucial to include the off diagonal elements in 
order to get the right stationary solution as we show in 
the following example. 

We adjust the voltages such that only transitions are 
possible from the ground state with 4m particles and en- 
ergy E\ m to ground and first excited states with 4m + 1 
particles and energies E\ mJrX and E\ m+l = El m+1 +eo, re- 
spectively. From the 4m + 1 ground states only transitions 
to the 4m ground state shall be allowed as depicted in Fig. 
[11] a). Note that this situation is only stable, because we 
have made the choice E c > eq. The master equation ex- 
pressed in the \n) basis reveals that only four of the 16 
states with the lowest particle hole excitation are indeed 
occupied in the stationary limit (cf. side b) of Fig. QT] 
because not all of the corresponding energetically allowed 
transitions from the TV = 4m ground state can be me- 
diated by one-electron tunnelling processes. Whereas in 
the |7V,m) basis, all 16 states with the energetically low- 
est bosonic excitations are equally populated. Since the 
degenerate states of the two bases are connected by a uni- 
tary transformation, the same must be true for the cor- 

I E 1 

responding matrix representations of Pq 4m+1 . From the 
representation of the RDM in the |n) basis, we know that 

IE 1 

the rank of p Q ' 4m+1 must be equal to 4. Because an uni- 
tary transformation does not change the rank of a matrix, 
the stationary solution in the | N, m) basis can maximally 
have 4 linearly independent columns. Since all diagonal el- 
ements are nonvanishing this is only possible if there are 
also nonvanishing coherences. 

Switching on the electron - electron interactions, the 
|n) states are no longer an eigenbasis of the SWNT Hamil- 
tonian and hence we must work in the \N,m) basis. But 
then, as we know from the discussion above, the coher- 
ences are expected to be of importance. 
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4.4.2 Negative differential conductance 

Spin charge separation and therefore non-Fermi liquid be- 
havior could also manifest itself in the occurrence of nega- 
tive differential conductance (NDC) at certain excitation 
lines involving transitions to states with fermionic exci- 
tations, as was predicted for a spinful Luttinger liquid 
quantum dot [24J with asymmetric contacts. We also find 
this effect for the nonequilibrium treatment of the SWNT 
quantum dot. Since in our case only the energy spectrum 
of the c+ mode depends on the interaction and the other 
three modes have the same energies as the neutral system, 
rather large asymmetries are needed in order to observe 
NDC. In Fig. [T2] we show the current across the first ex- 
citation line for transitions from N = 4m + 1 to N = 4m 
in the A = case. The corresponding trace in the Hb-Vb 
plane is indicated in the inset of Fig.[l2]a). Here the origin 
of the NDC is that some states with fermionic excitations 
have lower transition rates than nonexcited states, since 
due to the increased energy of the c+ modes less channels 
are available for transport. In Fig. [12] b) we show some 
of the excited states with N = 4m electrons which are 
responsible for the NDC, because their transition ampli- 
tudes to states with N = 4m + 1 electrons are reduced 
compared with the one of the 4m ground state. Apart 
from the asymmetries all other parameters are chosen as 
for Fig. [9] Only for asymmetries a — "fd/js larger than 
around 45, clear NDC features are seen. 



5 Conclusions 

In this article we have analyzed the linear and nonlin- 
ear current as a function of the gate and bias voltage 
across metallic SWNT quantum dots. The properties of 
the SWNT itself were derived from a microscopic model 
including electron - electron interactions. Exchange and 
related effects, which become relevant for small diameter 
SWNTs, have not been taken into account. The energy 
spectrum of a metallic SWNT, which can be obtained 
with the help of bosonization, turned out to be highly de- 
generate as a consequence of both fermionic and bosonic 
excitations. In the linear bias regime, the degeneracy of 
the groundstates leads to a characteristic pattern of con- 
ductance peaks depending on whether the two branches 
of the dispersion relation are aligned or not. Leaving the 
linear regime, asymmetry effects become relevant. Thus 
measurements of the current at low bias voltages, in the 
linear and nonlinear regime, in principle allow the sepa- 
rate determination of the source and drain tunnelling re- 
sistances as a function of the gate voltage. At higher bias 
voltages also excited states become relevant. The correct 
calculation of the nonequilibrium dynamics of the system 
then requires the inclusion of coherences in the reduced 
density matrix between degenerate states with bosonic ex- 
citations. At intermediate bias voltages there is a consid- 
erable deviation between the transport calculations with 
and without coherences. We emphasize that for a non- 
interacting system with unpolarized leads, coherences do 
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Fig. 12. a) Differential conductance as a function of the bias 
voltage for different asymmetries a = 7d/7s of the coupling to 
the leads. The trace in the fib-Vb plane across the boundary 
between regions A and B is indicated in the inset. Only for 
asymmetries larger than around 45 negative differential con- 
ductance occurs. All parameters here are chosen as for Fig. [9] 
except for fcsT a value of 0.026 meV was used, b) The tran- 
sitions responsible for NDC. i) In region A of the inset from 
a) only transitions between ground states are possible at low 
enough temperature. The transition rate from the N = 4m 
ground state to the N = 4m + 1 ground states is given by 47 s 
(see equation lf74|l ). ii) In region B additionally excited states 
become occupied. For some of the states with N = 4m elec- 
trons and fermionic excitations the transition rates to states 
with N = 4m + 1 electrons and neutral bosonic excitations is 
decreased compared to the ground state rate as a consequence 
of the larger energies of c+ excitations, which are not available 
for transport yet. Notice that we only show the most important 
types of transitions from A^ = 4m to A^ = 4m + 1 that take 
place in region B. Since we are considering a nonequilibrium 
situation other types of transitions are possible in principle. 



not have to be considered if expressing the reduced den- 
sity matrix in terms of Slater determinants, formed by the 
one electron wave functions of the noninteracting system. 
Another consequence of the electron correlations is the 
formation of a non-oscillatory spatial dependence of the 
tunnelling amplitudes along the nanotube axis. For tran- 
sitions between states with energetically low excitations 
we find a strong suppression of the tunnelling amplitudes 
near the SWNT ends. Furthermore we have addressed the 
influence of the tunnelling contacts on the transport. We 
have shown that extended contacts described as 3D Fermi 
gas do not lead to a polarization of the contacts with re- 
spect to the two branches of the dispersion relation. We 
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think that a further investigation of this point for other 
types of contacts is worthwhile. 

Useful discussions with S. Sapmaz and support by the 
DFG under the program GRK 638 are acknowledged. 



A Description of the tunnelling contacts 

In this appendix we discuss the dependence of the rates 
r, Ztmn)f M ' on t ne properties of the tunnelling contacts. 
In Section |4] we have introduced the parameters $i rr ' (e) 
with equation j58j) ■ Here we show the calculation of these 
parameters assuming a 3D electron gas in the leads. Fur- 
thermore, we require that the tunnelling region extends 
over several sites of the SWNT lattice. Starting from l|20p 
together with (fl6f and (1331) we can characterize the part 

of fe/j^f"" 1 " 1 that depends on the tunnelling contact I 
by the following expression, 

T lM « m (e) :=^sgn(FF') f d 3 r / 'dV2?(r)Tj(r')x 
ppi J J 

^2,4>* q {r)<Piq{r')Y^ <fs S n(F)fF(r)<pl gn(F , )f> F ,(r')x 

^fF.{x)) nk (4 IF , a {x')) . (79) 

The wave functions in the leads are denoted 4>i q and the 
sum 53 I extends over all q values that correspond to 
the energy e. At low enough energies, oscillations of the 
product {ipfFo{x)) nk (ipl, FI(7 (y) ) can be ignored along 

\ / k'm 

the length of the tunnelling interfaces. Using 1(54]) we thus 
can rewrite 1T791) as 



%,nkk' m {e) = L^2sgn(FF')vi FF/ (A)x 

FF' 

d 3 r J d s r%*(r)T l (r')A(x)A(x')J2K( r )^( r ') x 

J2 ^gn(F)f F (r)<p* sgn(F , )f , F ,{r') {4>fal) nk (fr><rl) k , m • 
rr 1 a 

(80) 

Here the factor v lFF >(A) = e ^i. d ^s g n(F-F')A takeg into 
account a possible phase shift due to the band mismatch 
A and we have used the notation ■f/jfo-z = ipfK o-(xi)/A(xi), 
where xi = 0, L for I — s, d. Now the parameters <Piff>{s) 
from 1(58]) are recovered by the relation 



%,nkk' m {e) = V $lrr> (e) V] {i>f<yl) nk ( V'fVi J , , 
— — \ / k'r 



In the next step we exploit that the Bloch waves V 5 sgn(_F)rF 
from equation ([28)1 are only nonvanishing around the po- 
sitions of the carbon atoms in the SWNT lattice. On the 
length scale of the extension of the p z orbitals all other 



quantities in <I>i? f < are slowly varying. Hence we can rewrite 
the integrals over x and y as a sum over the positions of 
the carbon atoms 

*iw(e) = C-^^sgn( J F 1 F')^F'(^)x 

L FF' 



^fpagn(F)rFfp> sga {F')r'F' X! T l* i x R,p) T l( x R' ,P') X 
p:p' R,R' 



4>lq{xR,p)<i>lq(xR' ,p') 



e iFR e- iFR \ (81) 



where we have defined xr p := R + t p . The constant C 
results from the integration over the p z orbitals. For the 
3D electron gas the wave functions 4>iq{ x ) are simply given 
by plane waves, 



<f>lq{x) 



a tqx 



Therefore we can easily perform the sum over the wave 
numbers q associated with the energy e, 



^4>iq{x)<t>lq{y) ~ 

dip / dcos9e iq ^ x -y\ cose 
o J -i 1\e\x-y\ 



47rsin(g| £ \x - y\) 



The previous expression is peaked around x = y. For 
q\ e > 1/ \x — y\ the correlations between x and y drop off 
fast. It means that if q\ e is larger than 1/ao, where ao is 
the nearest neighbor distance on the SWNT lattice, corre- 
lations between different carbon atom sites are suppressed 
such that we arrive at the following approximation: 



tf q (XR,p)4>lq(XR',p') Ri 47T(5H,R'(5p,p', 



which yields 



<Zw(e) = 4nC—x 



E 

FF' 



sgn{FF')v lFF , (A) V / psgll ( F ) rF / 



'psgn(F')r' F' X 



R 



m(x R . p )\ 2 e< F ~ F ') R A 2 (R x ). (82) 



Since we assume an extended tunnelling region, the fast 
oscillating terms with F = —F' are supposed to cancel. 
Furthermore we assume that both sublattices are equally 
well coupled to the contacts, such that the sum over R in 
([82| should approximately give the same result for p = 1 
and p = 2. Hence we can separate the sum over p from 
the rest. Since the Bloch waves ^ sgn (F)rF and ^ S gn(F)r' f 
are orthogonal to each other for r/r', the relation 

E fp s B n ( F ) r F fpsgn(F)r' F = $rr' 
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must hold, as it can be deduced from the explicit expres- 
sion for the Bloch wave, equation |28|) . and we arrive at 



${rf'(e) 



= <5 rr '47rC 



Nr 



\Ti(xR tP )\ A 2 (R X ). (83) 



R 



We should mention that the treatment of the contact ge- 
ometry here relies on several assumptions that might not 
be fulfilled under all circumstances and so it should be 
seen as a first estimate. For example we assume that the 
transmission functions TJ do not depend on the q vector of 
the incoming wave. But such a dependence would increase 
the correlation between different atom sites. Furthermore 
q\ e will not be much larger than l/a for all kind of con- 
tacts (for gold the Fermi wave number is about 1.5/do), 
nor can the lead electrons always be described by a 3D 
electron gas. 

Here we want to emphasize, the tunnelling contacts it- 
self play an important role for transport through a SWNT 
quantum dot. Especially it will be interesting for which 
kind of contacts the condition 

of having leads polarized with respect to the band degree 
of freedom f, might be fulfilled. 



where Xf^p(x) = - (a^. 9 f (x)) , since i<j)l <jF (x)+i(/)f r7F (x) 

is anti-hermitian. The actual values of the As are calcu- 
lated below in this appendix. Plugging the bosonization 
identity {84]) together with JSJ) into (N,m\ ip? aF (x) | AT' , m') , 
yields 



{N,m\^ f aF(x)\N',m') = 



K fF{N i) f<j {x)x 



( TO | e Z:,>o«- « /i, Ei,(45t(*)«<*«+^*(»)«J«,) | m ') . 

Remember that the factor (— 1)^j =1 Nj stems from the 
Klein factor. The term 

( TO | e i:,>o»- 0,/ "i: J «(^(»)«i*.+^'(*)«J,«) \ m ') (88) 

does not depend on the fermionic configuration and there- 
fore we have dropped the N index in \N,m). Now we 
exploit that in our finite size SWNTs there will always be 
only a finite number of bosonic excitations. Hence, by us- 
ing the Baker-Hausdorff formula, e A+B — e A e B e~^ A ' B ^ 2 , 
we commute all annihilation operators aj aq to the right 
and all creation operators a^ Sq to the left in 



B The matrix elements of the electron 
operators 

In this appendix we calculate the expressions for the ma- 
trix elements (N,m\ iprvF(x) \N',m') using the bosoniza- 
tion identity (|52|) of the electron operators ipfaF{x), namely 



1praF(x) 



\/l - e- a T 



(84) 



As a reminder, r\ ra is the Klein factor reducing the electron 
number by one, the operator K? F jsf f(T is given by 



KfFAf Tcr (x) :-- 



1 



2L 



giJ-agnF(sgnf-A/js<r+5)a; 



(85) 



and yields a phase factor depending on the filling of the 
band (fa). The bosonic fields 4>fFa(x) are given in terms 
of the bosonic annihilation operators b aq , 



i(j)faF(x) 



IT 

q>0 



isgn(Fr)qx—aq/2 



®(TSgn(r)q ■ 



(86) 



The action of ip rFcr (x) on the states \N, m) is conveniently 
determined by rewriting the operators b asgn ( r ) q in terms of 

the operators ajgq and Oj S which diagonalize the SWNT 
Hamiltonian Hq . As we know from the Bogoliubov trans- 
formation, the b operators depend linearly on the a oper- 
ators and hence we can write in general 

W? aF ( x ) + i<i>rc?F(x) = 

= E ^ E + Agg?(*)a}„) , (87) 



q>0 



( m | e E,>a«-° ,/ "E i ,(AKi-W»i.«+Ag,T(»)»J t «) \m') = 

n n (m j5q \e^ s ^H^M^ \ m .. 5q) t m 

q>0 jS 

— 1/2 / t \ m i 6K - 

with \mjs R ) = (rrijsJ) { a jS K J 1°) • Expanding the 

exponentials e****'^' 1 *** in (|89| . all terms which are of 
higher order than m'jgq w iU vanish. Analogously, all terms 

\ffijSq/ A t 

in the expansion of e faF ( - x,a ^i of higher order than rrijSq 
can be ignored. Hence we get 



(m|e Affia V a \m') =: F(X, m, to') 



E 



(A®at 



E 

3=0 



(XaY 



(90) 



Here, in favor of readability all indices have been dropped. 
In l|90p only the terms with m — i = m' — j survive. For 
to' > to we conveniently express j in terms of i and get: 



F(\,m,m') = A" 



nl \ 
7n l^i ' 



(A® A) 1 



in.' i i\(i + to' — m)\ (m — i)\ 



and for ml < to we write i in terms of j and get 
F(X, m, m) = 

(X®\) j ml 



(A®)' 



/ lit ! ^— \ 
m! ^— J 



ml j\(j + to — to')! (to' — j)! 



(91) 



(92) 



18 



Leonhard Mayrhofer, Milena Grifoni: Linear and nonlinear transport across carbon nanotube quantum dots 



Defining m max / min := max / min(m, to'), we can summa- 
rize J91| and {921) as 

F(X, to, to') = 

= (e(m' - m)X m '- m + 6>(to - to') (-A*) m ~ m ') 

( _ |A|2 ) ^max! 



E ilfi. 



^max- ~h Wmax ^min)' (^min 



or in terms of the Laguerre polynomials [25J, 



(93) 



m min' rm m „-m mm /|\|2\ 



F(A,TO,TO') = ,/^^i/ 

(0(m' - m)X m '- m + <9(to - to') (-A*) m ~ m ') . (94) 

Similar expressions for F(X, to, to') were found by Kim et 
al. [26] when investigating charge plasmons in a spinless 
Luttinger liquid quantum dot. In the end we obtain for the 
matrix elements {ipraF{x)) nm the following expression: 

{N,m\4, fGF {x)\N\m') = 

SN+e^N'i-l) 12 ^ N} K? F{N ')^(x)A(x)X 

IT II F ( A ^F ( X ) > m jSq > TO j5g) . 



9>0 i« 



where we have defined 



A(x) 



C The parameters A 

We still have to determine the values of the parameters 
X^plx) for our case. Using equation (|86|) we can express 

the sum + i(j> r<yF (x) in terms of the bosonic op- 

erators bg-g, 

i4>l aF (x) +i(j) raF {x) = 2^ — 

q >o V n i 

( _ p -isgn(Fr)qxi t i isgn(Fr)qxi \ /qc\ 

^ C °crsgn(r) 9 + 6 "asgn{r)qj , \ JJ / 

which in turn are related to the operators ajSq and aj 5? 
via equation (|45| . Hence 



^^(x) + i<j) ruF (x) 



E^E 



j<5 <?>0 V'"'.' 

— e 



^(^(S^aj^ + C^a,)] (96) 



By comparing (|96|) to ([871 the values of the As can be read 
off: 



With the coefficients Cj5 9 , Sj5 g and Ai 5 f from @5), (07 
and i|48|) we get for j<5 = c— , s+, s— , 



D isgi\{Fr)qx 



-A J6 



(97) 



and for the c+ mode we have 



i 



2./W a 




'~ q cos(qx) + i. /— sgn(Ff) sin(qx) I . (98) 
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